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ABSTRACT 


Basic heat transfer empirical and analytic equations are 
applied to a double envelope airship concept which uses 
heated Helium in the inner envelope to augment and/or 
control gross lift. The convective and conductive terms 
lead to a linear system of five eguations for the concept 
airship, with the nonlinear radiation terms included by an 
iteritive solution process. 

A FORTRAN program is used to perform the tedious 
calculations, and graphed results are obtained for the 
variables of interest. These indicate that a simple use of 
airship engine exhaust heat will give more than a 30% 
increase in gross airship lift. Possibly more than 100% 
increase can be achieved if a "steam injection" heating 


system, with associated design problems, is used. 
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I. NATURE OF THE PROBLEM 


A. CONTROL OT CITT 

In a fixed or rotary wing aircraft, the control 7٣ 
is normally a simple problem of changing engine power and/or 
angle of attack on the wing or rotor, and desired lift is 
quickly attained up to the aircraft's capability limits. 
Lighter than air ships, however, have a more inherent 
problem in increasing or decreasing lift. Typically, 
changes are produced by adding or valving lifting gas, 
ballasting by some heavy and readily obtained material such 
as sand or water, or employing engine exhaust based water 
recovery systems to compensate for fuel burnoff. 
Disadvantages of these are well known: availability and 
cost of lifting gas; ballasting material is not always 
available or easily transferred, and is a one way operation 
once airborne; and water recovery systems are bulky and 
heavy. 

Control of the lifting gas temperature could provide a 
more flexible and complete means of lift control. In 
addition, a significant increase in lifting gas temperature 


will ideally provide a proportionate increase in gross lift. 


B: HISTORY 
Using heated lifting gas to control or augment the Gr 
lift of an airship is an old idea, and an obvious extension 


T2 


of early hot air ballooning. Burgess, in his book "Airship 
Design" , lists this idea in his Airship Fallacies chapter 
due to his estimation on the rate of heat loss. He does, 
however, otíer “...gas cell construction in which convective 
losses were greatly reduced might possibly result ina 
practical solution of the problem" (reference 1, pp. 288- 
289). Present day fabrics, films and laminates offer much 
more flexibility in such construction than the rubberized 


cotton in use in the early 1900's (references 1 and 2). 


C. RESEARCH FINDINGS 

A library and information search was conducted to 
determine whether research results were available for such a 
heat transfer problem. The following results were obtained: 


- Mention in NASA CR137692 (reference 3) that Airfloat 
Transport, Ltd., had looked at heated Helium at takeoff 
to eliminate water recovery apparatus. No further 
information was obtained. 


- Concepts on using engine heat to heat lifting gas, but 
no technical evaluation by Davenport (reference 4). 


- A letter from Mr. I. P. Alexander of Wren Skyships, 
Ltd., stated that Wren Skyships had looked at ballonet 
heating on the RS-1. Additional mention was made of 
electrical heating tested on the LZ-130 Graf Zeppelin 
II. Neither of these concepts proved feasible, but 
these involved a small modification, and not a basic 
design concept. No technical data was received. 


- In a telephone conversation with Mr. Mark West of 
Aerostar, he stated that Aerostar had worked with 
students at Rensselaer Polytechnic Institute on a steam 
lift blimp design. No technical results were 
obtainable. 


- Extensive hot air and high altitude balloon heat 


transfer modelling is available. Applicable data from 
Stefan (reference 5) was modified for use in this model. 
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Also, applicable gas and film radiative properties were 
investigated in reference 6. 


The conclusion based on this available information is 
that only limited applications of lifting gas heating have 
been seriously explored. If a design dedicated to the use of 
heat to control and augment the lift of a Helium airship has 
been completed, no comprehensive mathematical modelling for 


such a design is available. 
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E THE PHYSICAL _ MODEL 


A. A PRACTICAL MODEL 

A comprehensive design concept for an airship is 
certainly beyond the scope of this paper. It is, however, 
desirable to have a reasonable physical concept on which to 
base the mathematical modelling and derivations. That 
concept, and some promising alternatives, are briefly 
explored, and the basis for the mathematical model is set 
Eorth. 

As encountered in any design, conflicting requirements 
must be considered, and the final choice of the model is an 
optimized compromise. Reduction of convective heat 
transfer, and maximum resistance to heat loss requires a 
multiplicity of panels within a volume, and insulative 
material coated onto panels. It is quickly realized, 
however, that such an airship will not be "lighter than 
air", and a minimum of envelope yardage and insulation must 
be pursued. 

In the absence of design specifications and "tradeoff" 
parameters, and in order to generate a simple but effective 
hypothetical physical model, it is assumed that the weight 
penalty of insulation material in addition to the inherent 
insulation properties of the envelope fabrics is 


prohibitive. Likewise, the use of additional fabric panels, 
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other than a basic double envelope design, is assumed 
prohibitive. The resulting concept is basically a 
conventionally shaped airship (see Figures 1 and 2) with the 
outer envelope enclosing lifting gas and ballonet volumes, 
providing basic protection from the external environment, 
aerodynamic shape, and pressure containment to hold that 
shape. 

An inner envelope contains the lifting gas, and must be 
expandable to accommodate both heating and altitude 
variations required for the airship. This envelope, then, 
must also carry the weight of the airship and loads. The 
weight suspension system also needs to constrain the shape 
of this inner envelope in such a way that the gap between 
inner and outer envelopes is reasonably uniform in order to 
minimize heat loss. The concept depicted in Figure 2 is of 
an adjustable/variable catenary cable system suspending the 
airship loads from the top of the inner envelope. The 
lengths of the cables are varied through a winch and pulley 
system to hold a desired pear shape, or natural buoyancy 
form, within the outer envelope. 

The air space between the envelopes is used as the 
ballonet in a standard fashion, and could, of course, be 
subdivided as required for safety and airship control. This 
volume also proves to be the most effective heat loss 


barrier as might have been expected. 
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Figure 2: Inner Envelope Suspension 
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In order to get a solid "feel" of the results from the 
heat loss rate calculations, specific external dimensions 
and gross lift capability were chosen to coincide closely 
with the probable sizes of near future NAVY or Coast Guard 
airships. These are: 

- gross lift of 130,000 lbf (2 x 10° ft3 of He) 
e length or 398.9 ft 

- maximum diameter 133 ft 

- slenderness ratio (L/2Rmax) = 3 

- total volume = 4 x 10° ft? 


Note that the outer envelope volume is twice the initial 


inner envelope volume. 


B. HEATING CONSIDERATIONS 

A practical, high capacity, efficient and lightweight 
system to provide heat into the inner envelope volume is 
critical to this concept. Low heat intensive systems such 
as heat pumps are assumed to be prohibitively heavy in the 
capacity ranges required. A simple air fed and vented 
petroleum based liquid burning heater may be a reasonable 
solution to these requirements. 

Davenport, in reference 4, proposes the use of engine 
exhaust to heat this volume. Appendix C shows a simple 
derivation of the amount of waste heat produced in engines 
such as may be used for propulsion of an airship this size. 
The problem of efficiently transferring this heat to the 
inner volume still remains. Forced air ducting and/or 


exhaust ducting may be the best, if least efficient means. 
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Since the heat exchange system required between hot 
ducting or heater surfaces and the inner Helium volume might 
be extensive and heavy, the injection of steam produced from 
a low pressure boiler or engine manifolds directly into the 
Helium is considered as a possible alternative. Section 
III,B,4 of this paper contains some comments on the 
difficulties of estimating the inner envelope heat transfer 
characteristics with saturated Helium and the condensation 
of water vapor on the inner envelope. Some foreseen 
detrimental factors in the use of such a system are the 
weight of the condensate on the inside of the inner 
envelope, corrosion or environmental effects on the inner 
envelope, weight of water required, and the energy to 
evaporate that water. A rough calculation of these values 
is completed in Appendix D. Benefits of the use of such a 


system are seen in section IV. 


C. ALTERNATIVES AND IMPROVEMENTS 

Using the inner volume as ballonet, with the lifting gas 
volume between the envelopes, has two significant advantages 
over the chosen model. First, the more elaborate suspension 
system for the inner envelope is not required. Second, an 
open burner or direct injection and mixing of heated air 
into the inner volume makes a heat exchange ducting network 
unnecessary. A third benefit, heating the heavier air, and 
thus reducing the density proportionately, is seen to be a 


fallacy due to the fact that a fixed Helium mass is 
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retained, while the air mass is discharged to hold a 
constant total volume. This alternative concept has the 
disadvantage of a reduced heated volume as temperature or 
altitude are increased, while in the chosen concept, the 
directly heated volume will increase with these parameters. 

An easy improvement to the modelled concept would be the 
use of lightweight panels such as a mylar film to break up 
the convective patterns. No estimate of the effectiveness 
of these panels is attempted in the present model. 

Further research might result in fabrication techniques 
for envelopes with dramatically increased heat loss 
resistance. One approach is to trap relatively thick air 
layers between mylar films of the fabric composites, which 


could result in a factor of 10 increase in this value. 


2 


III. THE MATHEMATICAL MODEL 


A. THE HEAT TRANSFER PROBLEM IN GENERAL 

Analytical solutions to the convective heat transfer 
problem are limited to laminar flow over specified 
geometries (reference 7). Empirical data curve fits must be 
used to apply known heat transfer characteristics to this 
airship heat loss problem. Since the available empirical 
equations apply to specific geometries, the airship shape is 
initially idealized as in Figure 3, with a hemispherical 
nose, cylindrical body, and conical tail. 

With a fixed slenderness ratio for the external shape, 
L/2Rmax = 3, the relationship between maximum diameter, 
envelope surface area and volume in this idealized model is 
fixed from basic solid geometry as: 

Rnax = (37/137)? = (4/7 e 
A=  (8+45)(3v/137 ) 2/3 

The inner envelope is idealized as a shape concentric to 
the outer envelope (Figure 3), with the total length the 
same as the outer envelope, but with the diameter varied as 
required to contain the lifting gas volume with changes in 
temperature (and altitude). With a constant proportion for 
the hemispherical nose and conical tail, the required length 
is derived for the cylindrical body, and again from solid 


geometry: 


22 


ano Z 





> 


5012115 5300134113 232571 


6 14 


<— noq 





23 


V ہر‎ ٠۷٦٤۵ 
Ain = 2ٔ ٣۴ = 3Ri n} 

In these idealized geometries, volume and area are 
greater than typical airship envelopes for a given maximum 
diameter. The more important volume and area relationship, 
however, is very close. These relationships are used 
because of their simplicity, allowing quick computer 
iterations, and the variation in manufactured airship 
envelope shapes and their corresponding geometric 


relationships. 


B. CONDUCTION AND CONVECTION MODELLING 
l. Conduction Through the Envelopes 

Since the dimensions of the airship are very large 
compared to the thickness (t) of the envelopes, a one 
dimensional heat conduction equation is completely valid 
(reference 8, ch. 1): 

q = -(KA/t) (Tout © Tin) 

The envelope thickness is idealized as a homogeneous 
material, although Vadala, in reference 2, clearly expects a 
multilayered fabric composite, such as a Kevlar base with 
covering Mylar films. The thermal conductivity (k) values 
used in this model are derived from reference 9 based on 
properties of Kevlar fabric in an appropriate thickness of 


.05 in. 


k = 0.03266 BTU ft/hr ft° F 
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2. External Convection 

Depending on operating conditions, the external flow 
will meet the criteria for forced, mixed, or free convection 
heat transfer conditions. In this approach, only forced and 
free convection values are derived, and a curve is sketched 
between zero and ten knot forced velocity values to 
represent the mixed convection regime. 

Operationally, this regime is expected to prove very 
important, as heat up and cool down periods will be required 
on takeoff and landing, depending on airship loading. With 
the airship tethered, light winds or even eddy currents in 
relatively calm air, with the boundary layer effects of the 
earth's surface, are assumed to make a direct analysis of 
this regime inaccurate. Therefore, even though a limited 
analysis of this mixed flow regime is possible, only the 
rough curvefit is used in the present model. 

For the free convection condition, the airship is 
idealized as a horizontal cylinder of equivalent area, and 
equation 7-36 from reference 8 is chosen as the closest 
empirical fit to the physical conditions: 


(0.6+0.387[GrqPr/[1+(.559/Pr)9/19]16/9]1/6,2 
where: k,«4, P, Pr, and Gr are based on the fluid 


temperature, 
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and Gr, is based on the maximum diameter "er the airship aces 
the temperature profile, Figure 4). 

Under forced convection conditions, the airship 
external surface is idealized as a simple flat plate with 
the equivalent length of the airship and the area of the 
outer envelope. Equations 5-79 and 5-83 of reference 8 are 
manipulated to yield the average equivalent convective heat 
transfer coefficient: 

h. = (K/L)Re;Pr!/3{[.2275/(log Re])%:584) 
- [850/Rer]) 
assumptions are: 


- laminar-turbulent boundary layer transition occurs at 
Re= 5x10 


- boundary layer is attached over the entire surface 


- the additional three dimensional body of rotation 
effects at the nose and tail are neglected 


- isothermal surface 
The derived convective heat transfer coefficients 
are then substituted into the linear heat transfer equation: 
qe hae OR 
3. Free Convection Between the Envelopes 
Here, due to the limitations of empirical data, the 
rough idealization of concentric horizontal cylinders is 
used. Equations 7-60 through 7-63, again from reference 8, 
are the closest approximation to the physical model 
available. These equations are based on isothermal 


Surfaces, and a constant spacing in the annulus between the 
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cylinders. A major problem here is that this empirical data 
is based on GrPr factors up to only 10Š. The airship model 
GrPr combination is on the order of 1013. However, since 
GrPr = 108 is well into a turbulent convective boundary 
layer, it is assumed there are no significant changes up to 
GrPr = 101}. Then: 

Kk, = K(0:4) [COB] os 


where: Gr, is based on the spacing between the cylinders, 


Rout ~ Rin = 4 


and the properties A, @ ,k are evaluated at T, 


Tn = the "volume mean" temperature 


3000000۴ ا ات رت رد تق )٢‏ 

Rm = (Rin + Rout)/2 
then: q= 2,7 Lko(T3 ات‎ 4111 Er 

4. Free Convection Inside the Inner Envelope 

Modelling the convection in the inner envelope 

proves to be the most elusive problem. Empirical eguations 
based on heat transfer from a uniform temperature fluid to a 
container wall are not adeguate due to the additional 
convective turbulence in this volume gained from whatever 
heat source is used.  Idealizing a heating duct and the 
envelope as concentric bodies results in Grashof numbers 
several orders of magnitude higher even than that between 
the envelopes. It is believed that completely different 


convective patterns will result, and that the application of 
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concentric body equations to this physical model is 
questionably accurate. 

The most applicable heat transfer coefficient 
available for this model can be derived from Stefan's 
article on hot air balloons (reference 5). Although the 
convective activity for the airship heater or heat source is 
less violent than the open flame of the hot air balloon, the 
scale of convective activity is of the same order, and the 
heat transfer mechanisms at the envelope wall are assumed to 
be similar. Accordingly, Stefan's convective heat transfer 
coefficient, based on hot air balloon flight data, is 
adapted for the Helium in the airship by assuming 

los موجہ‎ GrP 
as in the standard empirical forms. The exponent, n, is 
approximated as 0.25, similar to the value used in the 
volume between the envelopes. Then the constant is 
cancelled by: 
Nuys = Nusir (GEBE rs? / (GrPr) 41,277? 
and 
o b 


e eoa 


5 ( (Hair £He/MHe Cair)2]0:-23 /k. ; 


alr 
From reference 4, 
" 2 1 123 
DO St BTU/hu fe RIT — T2] / 


and 


MS... ..... 
This approach results in values of h double that 


derived from the application of the standard empirical 
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equations previously referenced. It is believed to be most 
accurate due to the similarity of convective heating action 
and size scale. 

Additional complications arise in this internal 
volume if a steam injection method is used as a heating 
source, as described in section II,B of this paper. Basic 
heat transfer texts such as references 8 and 10 discuss this 
difficult problem at length. Convective heat transfer will 
decrease due to the reduced temperature of wet steam from 
that of a simple liquid fuel heater, and a small insulative 
effect of the condensate film on the envelope is gained. 

The significant change, however, is seen when the heat of 
vaporization is transferred almost directly to the envelope 
wall during condensation. A direct estimate of the heat 
transfer characteristics is not attempted for this case, but 
data based on opposite extremes is obtained in section IV. 
The optimistic extreme is assumed as no increase in the heat 
transfer coefficient. The pessimistic extreme is the case 
where the envelope wall temperature is the same as mean 
fluid temperature (or an infinite heat transfer 


coefficient). 


C. RADIATION HEAT TRANSFER 
1. General 
Initially, the low temperature expected in this 
airship would lead one to neglect the radiation heat 
transfer. The insulative effect of the inter-envelope air 


DIO 


space, however, reduces the convective heat loss to the 
order of the radiative heat loss between the envelopes for 
expected envelope emissivity values. 

Envelope emissivity, then, becomes of primary 
importance. Developments in aerospace metallic coatings 
used to control temperature of satellites provides the 
promise of desired characteristics. Reference 11 discusses 
vacuum metallizing or solution processing techniques of 
depositing a thin metal coating onto mylar film which will 
result in an opague metal thickness of 5 x 1079 inches. 
This thickness results in about 25 pounds of Aluminum for 
each envelope surface coated, and emissivities less than 
o pale 

The addition of the non-linear radiative heat 
transfer term in a set of linear eguations becomes a major 
concern. Section D describes how the problem is handled in 
this model, and the problem of iteration convergence. 

2. Radiation to the Inner Envelope 

The straightforward solution of the radiative heat 
transfer between an enclosed body and the enclosing surface 
leads to equation 8-43 of reference 8: 

رر یا ےتک 
]1 — 10176 مر Mane La‏ 
where Abtr: Thtr. and Cina can be chosen appropriately. In‏ 
the results section, from program HOTSHIP, these values are‏ 


used pessimistically as: 
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A = SOE 

Thtr = 700 F 

Ente ~ 0-1 (typical of Puero 
and results are derived for a range of envelope 
emissivities. 

Computations are more complex for steam injection. 

Strong absorption lines for the water vapor lead to 
calculations of radiative heat transfer based on the mean 
free path computations as shown in chapter 8 of reference 8, 
and in reference 6. Due to the absence of the heater with a 
high surface temperature, however, the total radiative loss 
is expected to be less with the steam injection. This 
computation is not completed in this model, and program 
HOTSHIP assumes the radiative heat loss to be the same for 
both heating methods. 

3. Radiation Between the Envelopes 

The greatest radiative heat loss problem is between 

the envelopes, due to the large radiative surface areas. 
Again, the enclosed body equation is used: 

qr = DET gp 

(1/Ein)+[Ain/Ao][1/E5 - 1] 


and results are obtained for envelope emissivities from 0 to 





029% 


4. Radiative Loss from the Outer Envelope 
References 5 and 6 have a detailed treatment of 


radiation heat transfer between a hot balloon surface and 
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ene surrounding environment. A direct application to this 
airship model can be made, but this is not completed for 
three reasons. First, from the high convective heat loss 
only, the external envelope temperature is normally within 
two degrees F of the outside air temperature. Second, the 
expected operational conditions of this airship are assumed 
to vary to extremes, and results from reference 5 show that 
radiative heat loss is highly variable in these conditions. 
Accounting for these variations is not considered 
appropriate or necessary in this initial model. Third, due 
to expected night and foul weather operations, the external 
envelope emissivity should be small, minimizing the external 
radiative heat loss. Under specified environmental 
conditions, program HOTSHIP can be modified to include this 


external radiation term. 


D. THE RESULTING SET OF LINEARIZED EOUATIONS 

Discretizing the temperature profile as in Figure 4, 
then, the linearized set of heat transfer equations becomes 
a linear system of five equations in five unknowns for this 
particular physical model. 

At any surface, 

qd = Aconv * Ar 
In the inner envelope, then, 
do aoon FE 
3 ۹ 4 49 + ud ری‎ 
area, JA: (VE; —- 1] 
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and to linearize, an initial value is substituted for T, in 
the q, term, and Ir is treated as a constant in the linear 
system. Then, 

dg - hjAj(To => 12) + dro-2 


Similarly, between the envelopes, 


q= 2 7 kg L(T; - T,) + G Antr(Thtr 12) 
In (Rout/Rin) (1 htr)+[Ahtr/Ail[l ¡ 7 1] 
= 2 ZT s L(T4 - Tạ) + Ir3-4 
In (Rout/Rin) 


where initial values of T4 and dia ya substituted to make 
روب‎ a constant in the system. 


Rearranging terms, the resulting set of equations 


becomes: 
"iu q/h4A, = To + (Grg-2)/h3A, 
T>, - T} -q t;/KjA; = 0 
702 ۲ _ -q زا‎ E ال‎ 
e e 
T4 - Ts - g t,/kA, = 0 
T5 - م2۵‎ mm 
In matrix form, 
lo 0 ES 10> To + (9-02) 0008 
1-1 0 0 -t;/K¡A¡ JE 0 
0 1-1 0 -ln(Rout/Rin) Tu |= -d,-3-41n1(Rout/Rin) 
2 I Kol 2 W Kel | 
001.10... ۸۷ Ts 0 
oO O 1 A q Tinf 
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BE TIỀN 5: VALUES OBTAINED FOR INCREASE ENVELOPE HEAT LOSS 
RESISTANCE 

BOR A FACTOR OF 2 INCREASE: 

1 

EUT THE EMISSIVITIES OF THE INNER AND OUTER 

ENVELOPES, ELN AND EOUT 

? 


~ 


INPUT U INFINITY (UI,KTS), AND T INF (DEG R) 
: 

50,540 

NUMBER OF ITERATIONS ARE:RIN 5 ,FOR Q 8 
DODTN = 39.79 DMDT = 6.99 

RIN = 45.75 AIN = 102801.46 

VIN = 2115384.62 DVH20 = 254724 

PDL = 2.58 PDLTOT= 3.98 


ORE TEMPERATURE PROFILE (TO,T2,T3,T4,T5) IS: 
5500 548 . 3 548.47 540.41 540.16 

IMEMMETRIX VALUES A,B,C,D,E ARE: 

02032281 0.008935 0.553427 0.006463 0.003923 


DODTRO = 3.54 DỌDTR3 = 25.23 DQDIN= 39.79 
NUMBER OF ITERATIONS ARE:RIN 5 ,FORQ 9 

DODTN = 86.46 DMDT = 15.18 

RIN = 46.22 AIN = 103737.93 

VIN = 2153846.15 DVH20 = 79283.17 

PDL = 5.14 PDLTOT= 9.51. 


Moe TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS: 
560.00 557.82 557.06 540.90 540.34 

INMENINTIRIX VALUES A,B,C,D,E ARE: 

0.026212 0.008855 0.474857 0.006463 0.003923 


DODTRO < 3.40 DODTR3 = 52.43 DODIN = 86.46 
NUMBER OF ITERATIONS ARE:RIN 5 ,FORQ 9 

DODTN = 137.58 DMDT = 24.16 

RIN = 46.68 AIN = 104666.37 

VIN = 2192307.69 DVH20 = . 44 

PDL = 7.68 PDLTOT= 14 32 


ERESTEMPERATURE PROFILE (TO,T2,T3,T4,T5) IS: 
570.00 565:89 565.68 541.43 540.54 

Pee iat Rix VALUES A,B,C,D,E ARE: 

0.023150 0.008776 0.431857 0.006463 0.003923 


DODTRO = 3.24 DODTR3 = 81.42 DODTN = 137.58 
NUMBER OF ITERATIONS ARE:RIN 5 ,FORQ 10 | 
DODTN = 192.52 DMDT = 33.81 

RIN = 47.15 AIN = 105586.99 

VIN = 2230769.23 DVH20 = 192264.72 

PDL = 10.19 PDLTOT= 20.64 


HEM TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS: 
580.00 576.00 574232 542.00 540.76 
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THE MATRIX ۷5۲055 ۵ 8 0 5 SE 
0.020918 0.004350 0.399194 0.003231 0.003923 


DODTRO = 3.08 DQDTR3 = 116.51 DQDIN= 200.58 
NUMBER OF ITERATIONS ARE:RIN 5 ,FORQ 10 

DODTN = 262.07 DMDT = 46.02 

RIN = 47.61 AIN = 31 

VIN 7 = 2269230.77 DVH20 = 333175.90 2 
PDL = 12.65 PDLTOT= 30.69 


THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS: 
590700 584.95 583.82 541.87 541.03 

THE MATRIX VALUES A,B,C,D, E TAREE 

1 لات 2 0,00 0.376496 0.004312 0.019469 


DODTRO = 2.92 DODTR3 = 150.64 DODTN = 262.07 
NUMBER OF ITERATIONS ARE:RIN 5 ,FORQ 10 

DODTN = 327.36 DMDT = 57.49 

RIN = 48.07 AIN = 107405.63 

VIN = 2307692.31 DVH20 = 453004.31 

PDL = 15.10 PDLTOT= 39.47 


THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS: 
600.00 594.05 592.65 542.34 541.28 


THE MATRIX VALUES A,B,C,D,E ARE: 
0.018332 0.004276 0.358000 0.003231 0.003923 


DODTRO = 2.74 DODTR3 = 186.84 DODTN = 327.36 
NUMBER OF ITERATIONS ARE:RIN 5 „FOR O 10 

DODTN = 396.42 DMDT = 69.61 

RIN = 48.53 AIN = 108304.04 

VIN = 2346153.85 DVH20 = 655969.36 

PDL = 17.51 PDLTOT= 52.59 


THE TEMPERATURE PRORTLE (TO ,T2,.T2 4 +4٠٣٣٣ 
610,00 66071535 601.47 542.84 541.56 

IHE MATRIX VALUES A B C D E ARER 

0.017400 0.004241 0.342319 0.003231 0.003923 


DODTRO = 2.55 DODTR3 = 225.15 DODTN = 396.42 
NUMBER OF ITERATIONS ARE:RIN 5 ,FORQ 10 

DODTN = 469.23 DMDT = 82.40 

RIN = 48.98 AIN = 2 

VIN = 2384615.38 DVH20 = 1071527.20 

PDL = 19.89 PDLTOT= 76.87 


THE TEMPERATURE PROFILE (TO,T2,T3,T4,T5) IS: 
620.00 612.24 610.27 543.36 541.84 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.016613 0.004206 0.328650 0.003231 0.003923 
DODTRO = 2.36 DODTR3 = 265.63 DODTN = 469.23 
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VIN = 2423076.92 DVH20 = 1481890.78 

PDL = 21.91 PDEIOT= 100251 

THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS: 
630.00 389.08 31.00 571500 557.00 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.009450 0.004172 0.000000 0.003231 0.003923 


DODTRO = 2.84 DODTR3 = 0.00 DODTN = 4333.20 
NBMBER OF ITERATIONS ARE:RIN 4 ,FOR 0 9 

DODIN = 4897.97 DMDT = 860.13 

RIN = 49.89 AIN = 110957.35 

VIN = 2461538.46 DVH20 =  2272776.29 

PDL = 23.92 EDELOT- 144,22 


THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS: 
640.00 2 575.04 575.04 559.21 
THE MATRIX VALUES A,B,C,D,E ARE: 
0.009128 0.004139 0.000000 0.003231 0.003923 
DODTRO s 2.71 DODTR3 = 0.00 DODTN = 4897.97 


NUMBER OF ITERATIONS ARE:RIN 4 ,FOR Q 9 


DODTN = 5472.35 DMDT = 961.00 

RIN = 50.33 AIN = 111828.94 
VIN = 2500000.00 DVH20 = 4419910.85 
PDL = 26.10 PDLTOT= 258081 


THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS: 
650.00 601.63 579.15 579.15 561.47 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.008844 0.004107 0.000000 0.003231 0.003923 


DODTRO = 2.59 DODTR3 = 0.00 DODTN = 5472.35 


NUMBER OF ITERATIONS ARE:RIN 5 „FOR O 9 


DODIN = 6055.73 DMDT = 1063.44 
RIN = 50.78 AIN = 112694.08 
VIN = 2538461.54 DVH20 = 8108853.38 
PDL = 28.36 PDLTOT= 452.45 


THE TEMPERATURE PROFILE (TO,T2,T3,T4,T5) IS: 
660.00 608.00 583.32 583.32 563.76 
THE MATRIX VALUES A,B,C,D,E ARE: 
0.008590 0.004075 0.000000 0.003231 0.003923 
DODTRO = 2.45 DODTR3 = 0.00 DODTN = 6055.73 


NUMBER OF ITERATIONS ARE:RIN 5 ,FORQ 9 
DODTN = 6647.57 DMDT = 1107.38 
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SECTION 3: VALUES OBTAINED FOR T3 = T4 


INPUT THE EMISSIVITIES OF THE INNER AND OUTER 


ENVELOPES, EIN AND EOUT 
? 


S 
INPUT U INFINITY (UI,KTS), AND T INF (DEG R) 
3 

50,540 

NUMBER OF ITERATIONS ARE:RIN 5 ,FOR Q 7 
DODTN = 332.86 DMDT = 58.45 

RIN = 45.75 AIN = 102801.46 

VIN = 2115384.62 DVH20 = 25472.64 

PDL = 2.32 PDLTOT= 3.74 


THE TEMPERATURE PROFILE (TO,T2,T3,T4,T5) IS: 
550.00 7 542.38 542.28 21 
THE MATRIX VALUES A,B,C,D,E ARE: 
0.018624 0.004468 0.000000 0.003231 0.003923 
DODTRO = 3.62 DODTR3 = 0.00 DODIN = 332.86 


NUMBER OF ITERATIONS ARE:RIN 5 ‚FOR © 1 


DODIN = 151.03 "DDI = 131299 

RIN = 46.22 AIN = 103937253 
VIN = 2153846.15 DVH20 = 29۲ 
PDL = 4.73 PDLTOT= 2 


THE TEMPERATURE PROFILE (TO,T2, +4۰۲ 
560.00 548.71 545.38 545.38 542,95 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.015099 0.004427 0.000000 0.003231 0.003923 


DODTRO = 3.54 DODTR3 = 0.00 DODTN = 751.63 
NUMBER OF ITERATIONS ARE:RIN 5 „FOR 0 8 

DODTN = 1207.26 DMDT = 212 01 

RIN = 46.68 AIN = 104666.37 

VIN = 2192307.69 DVH20 = 121494.44 

PDL = 7.16 PDLTOT= 13.86 


THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS: 
570.00 553.93 548.64 ' 548.64 544.74 
THE MATRIX VALUES A,B,C,D,E ARE: 
0.013346 0.004388 0.000000 0.003231 0.003923 
DODTRO = 3.46 DODTR3 = 0.00 DỌDTN = 1207.26 


NUMBER OF ITERATIONS ARE :RIN 5 TOD 8 


DODTN = 1688.06 DMDT = 296.44 

RIN = 47.15 AIN = 105586.99 
VIN = 2230769.23 DVH20 = 192264.72 
PDL = 9.60 PDLTOT= 20.14 


THE TEMPERATURE PROFILE (TO,T2,T3,T4,T5) IS: 
580.00 559 2 552.05 552.08 546.62 
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RINE = 21.22 AIN = 113552.84 
Will “= 2576923.08 DVH20 = 35853817.74 
ELE = 35199 PDLTIOT= 1359532 


THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS: 
970.00 670.00 565 2 547 . 2 543 . 6 
COS MATRIX VALUES A,B,C/D,E ARE: 
) 10000 0.004045 0.272200 0.003231 0.003923 
DODTRO = 0 9۹2 DODIR 2 513.25 DODIN = 1009.32 


MEMBER OF ITERATIONS ARE:RIN 5 „FOR O 6 


DODTN = Irro S DUDE = 196263 

RIN = 51.66 AIN = 114405.53 
MeN 2515584162 DVA20 = xx=xx=x=xxx=x 
PDL = 34.20 PDLTOT= - 752289 


THE TEMPERATURE PROFILE (TO,T2,T3,T4,T5) IS: 
680.00 680.00 675.50 548.01 544.39 
THE MATRIX VALUES A,B,C,D,E ARE: 
0.000000 0.004014 0.264220 0.003231 0.003923 
DODTRO = 0.62 DODTR3 = 637.19 DODTN = 1119.72 
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NUMBER OF ITERATIONS ARE:RIN 5 „FOR 0 6 


DODTN = 616.34 DMDT = 108.24 
RIN = 49.44 AIN = 110079.94 
VIN = 2423076.92 DVH20 = 1481890.78 
Pp = 22.79 PDLTOT= 100.52 


THE TEMPERATURE PROFILE (ZO TZ TT TS TS 
630.00 630.00 627.43 544.41 542.42 

THE MATRIX VALUES AVE 6:5 ake. 

0.000000 0.004172 O. 3103422 07003231. 0.005525 


DODTRO = 1.96 DODTR3 = 348.83 DODTN = 616.34 
NUMBER OF ITERATIONS ARE:RIN 4 „FOR O 6 

DODTN = 707.53 DMDT = 124.25 

RIN = 49.89 AIN =  ټ‎ ,.85 

VIN = 2461538.46 DVH20 = 2272776.29 

PLE 25.15 PDLTOT- 143.52 


THE TEMPERATURE PROFILE (TO,T2,T3,T4,T5) IS: 
640.00 0 637.07 545.06 542.78 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.000000 0.004139 0.299614 0.003231 0.003923 


DODTRO = 1.71 DQDTR3 = 400.44 DODTN = 707.53 
NUMBER OF ITERATIONS ARE:RIN 4 ,FORQ 6 

DODTN = 803.35 DMDT = 141.08 

RIN = 50.33 AIN = 111828.94 

VIN = 2500000.00 DVH20 = 4419910.85 

PDL = 27.46 PDLTOT= 250209 


THE TEMPERATURE PROFILE (T0,T2,T3.74,T5) 75: 
650.00 650.00 646.70 545.75 543.15 

IHE MATRIX VALUES A B,CD E AREF 

0.000000 0.004107 0.259787 0.003231 0 32.. 


DODTRO = 1.46 DODTR3 = 454.98 DODTN = 803.35 
NUMBER OF ITERATIONS ARE:RIN 5 ,FORQ 6 

DODTN = 903.90 DMDT = 158.73 

RIN = 50.78 AIN = 112694.08 

VIN = 2538461.54 DVH20 = 8108853.38 

PDL = 29.74 PDLTOT= : 446.08 


THE TEMPERATURE PROFILE (TO,T2,T3,T4,T5) IS: 
660.00 660.00 656.32 546.47 543.55 
THE MATRIX VALUES A,B,C,D,E ARE: 
0.000000 0.004075 0.280690 0.003231 0.003923 
DODTRO = 1.19 DODTR3 = 512.55 DODTN = 903.90 


NUMBER OF ITERATIONS ABE +-+ 70 6 
DODIN = 1009.32 1101 - i zo 
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SECTION 2: VALUES OBTAINED FOR T2 = TO 


INPUT THE EMISSIVITIES OF THE INNER AND OUTER 
ENVELOPES, EIN AND EOUT 
? 


0+ 00۹9 


INPUT U INFINITY (UI,KTS), AND T INF (DEG R) 
? 

50,540 

NUMBER OF ITERATIONS ARE:RIN 5 „FOR O 5 
DODTN = 47.35 DMDT = 8.32 

RIN = 45.75 AIN = 102801.46 

VIN = 2115384.62 DVH20 = 25472.64 

PRO = 2.67 PDLTOT= 4.07 


MW TEMPERATURE PROFILE (TO,T2,T3,T4,T5) IS: 
550.00 550. 00 549.79 540.34 540.19 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.000000 0.004468 0.535970 0.003231 0.003923 


DODTRO = 3.52 DODTR3 = 29.72 DODIN = 47.35 
NUMBER OF ITERATIONS ARE:RIN 5 ,FORQ 5 

DODTN = 102.26 DMDT = 17.96 

RIN = 46.22 AIN = 3 

VIN = 2153846.15 DVH20 = 79283.17 

sep = 5.31 PDLTOT= 9.67 


BER TEMPERATURE PROFILE (TO ,T2,T3,T4,T5) IS: 
560.00 560.00 299.99 540.73 540.40 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.000000 0.004427 0.460701 0.003231 0.003923 


DODTRO = 3.36 DODTR3 = 61.42 DODIN = 102.26 
NUMBER OF ITERATIONS ARE:RIN 5 ,FORQ 5 

DODTN = 162.32 DMDT = 28.51 

RIN = 46.68 AIN = 104666.37 

VIN = 2192307.69 DVH20 = 121494.44 

PDL = 7.92 PDLTOT= 14.54 


RE TEMPERATURE PROFILE (TO,T2,T3,T4,T5) IS: 
570.00 570700 2 9 541.16 540.64 

ENE MATRIX VALUES A,B,C,D,E ARE: 

0.000000 0.004388 0.419335 0.003231 0.003923 


DODTRO = SEELE DOETR3 = - 8525 DODIN = 162.32 
NUMBER OF ITERATIONS ARE:RIN 5 ‚FOR 0 6 

DODTN = 227.01 DMDT = CE 

RIN = 47.15 AIN = (+0 s2 

Willie 2230769.23 DVHZO - £192264.72 

PDL = 10.49 PDLTOT= 20.90 


ERENTEMBERATURE PROFILE (T0,T2,T3,T4,T5) 15: 
580.00 580.00 5 001 541.62 540.89 
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THE MATRIX VALUES A,B,C,D,E ARE: 
0.013889 0.004045 0.277434 0.003231 0.003923 
DODTRO = 1.26 DODTR3 = 503.05 DODTN = 891.11 


NUMBER OF ITERATIONS ARE:RIN 5 „FOR O 12 


DODTN = 987.43 DMDT = 173.40 
RIN = 51.66 AIN = 114405.53 
VIN = 2615384.62 DVH20 = 大 大 大 六 六 大 大 大 大 大 大 
PDL = 33.48 PDLTOT= -809.11 


THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS: 
680.00 666.70 662.73 547.06 543.87 
THE MATRIX VALUES A,B,C,D,E ARE: 
0.013485 0.004014 0.269284 0.003231 0.003923 
DODTRO = 1.01 DODTR3 = 557.89 DODIN = 987.43 
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vapor at saturation level in the Helium becomes more and 
more rapid as temperature is increased. As previously 
noted, the loss of insulative space from the water vapor 
volume is not accounted for in program HOTSHIP, and the heat 
loss rates shown are increasingly optimistic with 
temperature rise. This effect is not more than a few 
percent at 120 deg F (delta T = 60) but by 160 deg F, the 
inner volume would be greater than the outer volume, 
implying that by this point the inter envelope insulative 
effect is absent, and the heat flux will actually be many 
times higher as shown previously in this section. 

If losses due to condensation on the inner envelope are 
much higher than expected, the pessimistic approximation of 
ene 15 made by setting the matrix value of "A" to zero. 
Results in Appendix B show only a 12.9% increase in heat 


Moss for this condition. 


B. ESTIMATES OF ACCURACY 

The best accuracy is obtained from the conduction 
equations in the envelope materials. For the material 
described, the conductivity is estimated to be within 20% of 
the average of that to be obtained in an actual envelope. 
The one dimensional conduction equation is completely 
applicable. Unfortunately, this is among the least critical 
conditions due to the ineffective insulation of the 
envelopes. Even if the conductivity of the envelopes was 
halved, the effect on heat loss is minor. In Appendix B, 
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Although some heat loss resistance is exhibited at the 
boundary layers on the inside and outside, the results file 
in Appendix B clearly shows that the inter-envelope volume 
provides most of the insulative effects. More than 8/10 of 
the temperature drop is seen across this air gap under 
"normal" cruise conditions, implying that the heat loss 
would be several times higher in a conventional airship 
design. In fact, if the matrix value "C" is set to zero in 
program HOTSHIP, the heat flux is seen to increase by a 
factor of 8.4 (see Appendix B). This result does not 
account for the additional radiative heat loss which will 
occur due to the significant temperature increase of the 
outer envelope. 

This insulative gap is seen to be particularly critical 
as alrship speed is increased from zero. Figure 10 shows 
that the heat flux quickly rises to an asymptotic value 
about twice that for zero velocity even in the concept 
design. This rise would be accentuated by the appropriate 
factor of more than eight without the inter envelope space. 
This implies that a standard airship moving at 50 knots will 
have an order of magnitude higher heat loss than a hot air 
balloon with the same shape and size, and 17 times the heat 
loss of the concept design at zero velocity. 

Figures 8 and 9 show a dramatic difference in lift 
increase (for the same heat input) when the steam injection 


heating concept is used. The volume increase due to water 
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400 


500 


naa 


600 


FORCED CONVECTION PARAMETERS OUTSIDE THE AIRSHIP, BASED ON A 
CRITICAL REYNOLDS NUMBER RECRIT=500000 
RHOI=0.0022846*540.0/TI 
MUI=((4.3339E-7-3.85605E-7)/90.0)*(TI-540.0)+3.85605E-7 
REL=RHOI*UI*6076.0*L/(MUI*3600.0) 
PRI=((.697-.708)/90.0)*(TI-540.0)+.708 
KINF=(((.0173513-.0151615)/90.0)*(TI-540.0)+.0151615)/3600.0 


H4=((0.2275/((LOG10(REL))**2.584))-850.0/REL)*KINF*REL* 


+(PRI**0.333333333)/L 


SETTING TRE INITIAL VALUES ۶۷ THE LINEAR SISIEN ٣٣ 


T2=T0-(T0-TI)/5.0 
T3=T0-(T0-TI)/4.5 
T4=T0-(T0-TI)*4.0/5.0 
T5=T0-(T0-TI)*4.5/5.0 
DODT=480 

EHTR=0.1 

AHTR=500 

THTR=700 


Jsl 


IF (UI.NE.0.0) GO TO 600 
IF EXTERNAL FLOW IS FREE CONVECTION, COMPUTE THE HEAT. TRANSIER 
COEFFECIENT FOR UI=0, IF NOT, R4 ]5 ھ۶٦‎ + + ٥ 


TFL=(T5+TI)/2 

PRFL=((.697-.708)/90.0)*(TFL-540.0)+.708 
KFL=(((.0173513-.0151615)/90.0)*(TFL-540.0)+.0151615)/3600.0 
MUFL=((4.3339E-7-3.85605E-7)/90.0)*(TFL-540.0)+3.85605E-7 
RHOFL=0.0022846*540.0/TFL 
GRFL=32.174*(T5-TI)*8.0*(ROUT**3)*(RHOFL**2)/(TFL*(MUFL**2)) 


H4=(KFL/L)*((.6+.387*(GRFL*PREL/(1.0+(.559/PRFL)**.5625)**1.7778) 


+**,1666667)**#2) 


FREE CONVECTION BETWEEN THE ENVELOPES BASED ON HORIZONTAL 
CYLINDRICAL ANNULUS, ISOTHERMAL ENVELOPES, AND CONSTANT 
ENVELOPE SPACING 


RMEAN=(ROUT+RIN)/2.0 
TMEAN=(( (RMEAN**3 ) - (RIN**3 ) ) *T3+( (ROUT**3 ) - (RMEAN**3 ) ) *T4)/ 


+((ROUT^^3)=(RIN^^3)) 


RHOMN=0 .0022846*540.0/TMEAN 
MUMN=3.85605E-7+((4.3339E-7-3.85605E-7)/90.0)*(TMEAN-540.0) 
PRMN=0.708+((.697-.708)/90.0)*(TMEAN-540.0) 
KMEAN=(0.0151615+((.0173513-.0151615)/90.0)*(TMEAN-540.0))/3600.0 
GRDz232.174*(T3-T4)* ( (ROUT-RIN) **3) * (RHOMN**2) / (TMEAN* (MUMN**2) ) 
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ME MAIRIS VALUES A,B,C€,D,E ARE: 
001 3/ 0.008700 0.402213 0.006463 3 


DODTRO < 3.08 DQDTR3 = 112.16 DODTN = 192.52 
NUMBER OF ITERATIONS ARE:RIN 5 ,FORQ 10 

DODTN = 251.11 DMDT = 44.10 

iN = 47.61 AIN = 106500.01 

VIN = 2269230.77 DVH20 = 333175.90 

PDL = 12.66 PDLTOT= 30.70 


THE TEMPERATURE PROFILE (TO,T2,T3,T4,T5) IS: 
390.00 5059.12 55 5 542.61 540.99 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.019680 0.008625 0.379454 0.006463 0.003923 


DODTRO = 2.91 DODTR3 = 144.80 DODTN = 251.11 
NUMBER OF ITERATIONS ARE:RIN 5 „FOR O 10 

DODTN = 313.19 DMDT = 55.00 

RIN = 48.07 AIN = 107405.63 

VIN = 2307692.31 DVH20 = 453004.31 

HE = 15.11 PDLTOT- 39.48 


THE TEMPERATURE PROFILE (TO,T2,T3,T4,T5) IS: 
600.00 594.24 591.57 543.25 541.23 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.018538 0.008552 0.360910 0.006463 3 





DODTRO = 2.74 DODTR3 = 179.32 DODIN= 313.19 
NUMBER OF ITERATIONS ARE:RIN 5 ,FORQ 10 

DODTN = 378.69 DMDT = 66.50 

RIN = 48.53 AIN = 4 

VIN = 2346153.85 DVH20 = 655969.36 

PDL = 17.52 PDLTOT= 52.59 


EM TEMPERATURE PROFILE (TO,I2,I13,14,15) IS: 
210700 603.38 600.17 543.93 541.49 

MUEMMATRIX"VALUES A,B,C,D,E ARE: 

0.017602 0.008481 0.345190 0.006463 0.003923 


DODTRO = 2.55 DODTR3 = 215.78 DODTN = 378.69 
NUMBER OF ITERATIONS ARE:RIN.5 ,FOR Q 11 

DODTN = 447.65 DMDT = 78.61 

RIN = 48.98 AIN = 109195.42 

VIN = 2384615.38 DVH20 = 1071527.20 

PDL = 19.91 PDLTOT= 76.87 


THE TEMPERATURE PROFILE (TO,T2,T3,T4,T5) IS: 
620.00 6402.51 608.75 544.65 541.76 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.016819 0.008412 0.331477 0.006463 0.003923 


DODTRO = 2.36 DODTR3 = 254.29 DODTN = 447.65 
NUMBER OF ITERATIONS ARE:RIN 5 ,FOR Q 11 

DODTN = 520.00 DMDT = 91.32 

RIN = 49.44 AIN = 110079.94 
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VIN = 2423076.92 DVH20 = 1481890.78 

DA 22.25 PDLIOT= 1O0 51 

THE TEMPERATURE PROFILE CIO TE TS aD IS. 
630.00 621.64 51.41 545.40 542.04 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.016136 0.008344 0.319288 0.006463 0.003923 


DODTRO = 2.15 DODTR3 = 294.80 DODTN = 520.00 
NUMBER OF ITERATIONS ARE:RIN 4 ,FOR Q 11 

DODTN = 595.80 DMDT = 104.63 

RIN = 49.89 AIN = 110957.85 

VIN = 2461538.46 DVH20 = 2272776.29 

PDL = 24.57 PDLIOT= 143.82 


THE TEMPERATURE FRORTER TO TI a) 
640.00 630.77 625.84 54638 542.34 

THE MATRIX VALUES A B,C, D,E ARES 

0.015535 0.008278 0.308269 0.006463 0.003923 


DODTRO = 1.94 DODTR3 = 337.41 DODIN = 595.80 
NUMBER OF ITERATIONS ARE:RIN 4 ,FOR Q 11 

DODTN = 675.08 DMDT = 118.55 

RIN = 50.33 AIN = 4ه‎ 

VIN = 2500000.00 DVH20 = 4419910.85 

PDL = 26.85 PDLTOT= 257.31 


THE TEMPERATURE PROFILE (10 12151... 7... 
650.00 639.90 634.35 547.01 542.65 
THE MATRIX VALUES A,B,C,D,E ARE: 


0.015000 0.008214 0.298185 0.006463 0.003923 
DODTRO = 1.72 DODTR3 = 382.16 DODTN = 675.08 


NUMBER OF ITERATIONS ARE:RIN 5 „FOR O 11 


DODTN = 757.89 DMDT 2 133500 

RIN = 50.78 AIN = 112694.08 
VIN = 2538461.54 DVH20 = 8108853.38 
PDL = 29,10 PDETOIS 449.04 


THE TEMPERATURE PROFILE (TO,T2,T3,T4,T5) IS: 
660.00 649.02 642.84 547.87 542.97 
THE MATRIX VALUES A,B,C,D E TARE E 
0.014518 0.008151 0.288857 0.006463 0.003923 
DODTRO = 1.49  DODTR3 z 429.12 DODIN = 757.89 


NUMBER OF ITERATIONS ARE:RIN 5S F CP : 1> 


DODTN = 844.26 DMDT = 148.26 

RIN = 51.22 AIN = 113552.84 
VIN = 2576923.08 DVH20 = 35853817.74 
PDL = 31.32 PDLTOT= 1875265 


THE TEMPERATURE PROFILE (10,22,73 727 pelo: 
670.00 65815 691,530 548.77 543.31 
THE MATRIX VALUES A,B,C,D,E ARE: 
0.014076 0.008089 0.280164 0.006463 0.003923 
DODTRO = 1.24 DODTR3 = 478.27 DODTN = 844.26 
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NUMBER OF ITERATIONS ARE:RIN 5 „FOR O 12 


DODTN = 934.32 DMDT = 164.08 
RIN = 51.66 AIN = 11440553 
VIN = 2615384.62 DVH20 = 大 大 大 大 大 大 六 大 大 大 大 
PDL = 33.50 PDLTOT= -808.84 


THE TEMPERATURE PROFILE (TO,T2,T3,T4,T5) IS: 

680.00 667.24 659.74 549.70 543.67 
IZEMMATRIX VALUES A,B,C,D,E ARE: 

512674 0.008029 0.271992 0.006463 “0.003923 
DODTRO = 0 95 bo Li = 529.77 DODTN = 934.32 


FOR A FACTOR OF 10 INCREASE: 


RMPUT THE EMISSIVITIES OF THE INNER AND OUTER 
ENVELOPES, EIN AND EOUT 


2 3 

NIN, .O 

INPUT U INFINITY (UI,KTS), AND T INF (DEG R) 
? 

50,540 

NUMBER OF ITERATIONS ARE:RIN 5 ,FOR Q 3 
DODTN = 31.55 DMDT = 5,54 

RIN = 45.75 AIN = 102801.46 

VIN = 2115384.62 DVH20 = 25472.64 

PDL = 2.60 PDLTOT= 4.00 


CEE TEMPERATURE PROFILE (TO,T2,T3,T4,T5) IS: 
290.00 549.03 547.62 541.14 540.12 

BUE MATRIX VALUES A,B,C,D,E ARE: 

0.034562 0.044676 0.577787 0.032314 0.003923 


DODTRO = 3.54 DODTR3 = 20.33 DODTN = 31.55 
NUMBER OF ITERATIONS ARE:RIN 5 ,FORQ و‎ 

DODTN = 67.25 DMDT = DT 

PIN = 46.22 AIN = 103737.93 

VIN = 2153846.15 DVH20 = 79283.17 

BDL = 5.17 PDLTOT= 9.54 


MEN TEMPERATURE PROFILE (TO,T2,T3,T4,T5) IS: 
560.00 5971 53 542.44 540.26 

IRE IINIRIX VALUES A,B,C,D,E ARE: 

0.028005 0.044273 0.497466 0.032314 0.003923 


DODTRO = 3.39 DODTR3 = 41.52 DODIN = 67.25 
NUMBER OF ITERATIONS ARE:RIN 5 ,FOR 0 10 

DODTN z 105.57 DMDT = 18.54 

RIN = 46.68 AIN = 7 

VIN = 2192307.69 DVH20 = 121494.44 

PDL = 7.72 PDLTOT= 14.36 


THE TEMPERATURE PROFILE (TO,T2,T3,T4,T5) IS: 
570.00 567.47 5642.03 543.83 540.41 
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NUMBER OF ITERATIONS ARE:RIN 53 12 


DODTN = 1013.73 DMDT = 178.02 

RIN = 51.66 AIN = 23 
VIN = 2615384.62 DVH2O = 大大 六 大大 大 大 大 大 六 大 
PDL = 33.21 POLIGI: -812.60 


THE TEMPERATURE PROFILE (TO,T2,T3,T4,T5) IS: 
680.00 666.43 662.36 543.28 540.00 
THE MATRIX VALUES A,B,C,D,E ARE: 
0.013397 0.004014 0.267670 0.003231 0.000000 
DODTRO = 1.02 DODTR3 = 568.82 DODIN = 1013.73 


23 


BEESMATRIX VALUES A,B,C,D,E ARE: 
0.024750 0.043880 0.453521 0.032314 0.003923 


DODTRO = 3223 DODIR -MES come DODIN = 105.57 
NUMBER OF ITERATIONS ARE:RIN 5 „FOR O 10 

DODTN = 146.21 DMDT = 25.68 

RIN = 27:15 214 = 9 

EI = 2230769.23 DVHZO = | 192264.72 

PDL = 10.24 PDLTOT= 2069 


DIES TEMPERATURE PROFILE (TO,T2,T3,T4,T5) IS: 
280.00 579.76 570.40 545 0 540.57 


— ey, 


Poe eMATRIX VALUES A,B,C,D,E ARE: 
0022670 0.043498 0.423124 0.032314 3 


DODTRO = 3.07 "008155 > 56.89 DODTN = 146.21 
DIDHBER OE ITERATIONS ARE:RIN 5 ,FOR Q 11 

DODTN = 188.96 DMDT = SS 

RIN = 47.61 AIN = 106500207 

DVHZO = 333175.90‏ 2269230.77 جح کت 

PDL = 12.73 PDLDTODS 30.75 


MEMTEMPERATURE PROFILÉ (T0,T2,T3,T4,T5) IS: 
290.00 8536.06 3 3 ٤۶ 546 5 540.74 

Mme MATRIX VALUES A,5B,C,D,E ARE: 

۲۶۰210164 0.043125 0.399790 0.032314 3 


DODTRO = 25958 0185 > 111:25 01111506 
NUMBER OF ITERATIONS ARE:RIN 5 ,FOR Q 11 

DODTN = 233.565 DIDIT = 41.03 

RIN = 48.07 AIN = 107405.63 

wane = 2307692.31 DVHZO = | 453004.31 

PDL = 175.18 PDETOT= 39253 


IH EMTENPERATURE PROFILE (TO,T2,T3,T4,T5) IS: 
600.00 595.39 585.40 548.47 540.92 

IESMMATRIX VALUES A,B,C,D,E ARE: 

02019974 0.042761 0.380820 0.032314 3 


DODTRO = 2271 DODTR3 2 136:58 DODIN =5 293.605 
MEMBER OF ITERATIONS ARE:RIN.S „FOR O 11 

DODTN = 280.24 DMDT = 49.21 

RIN = 48.53 AIN = 105304 . 4 

VIN = 2346153.85 DVH20 = | 655969.36 

PDL = 1760 EDITOIS 52.64 


DEESSEMPERATURE PROFILE (TO,T2,T3,T4,T5) IS: 
610.00 604.72 592.84 5505 541 0 
TRE MATRIX VALUES A,B,C,D,E ARE: 
0.019002 0.042406 0.364745 0.032314 0.003923 
DODTRO = 2.52 DODTR3 = 163.22 DODTN = 280.24 


NUMBER OF ITERATIONS ARE:RIN 5 ,FOR O 11 
DODTN 328.70 DHP 57012 
RIN 48.98 AIN 109195.42 
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VIN = 2384615.38 DVH20 = 1071527229 

EDI = 19.99 PDLTOT= 76 . 90 

THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS: 
620.00 614.07 600.24 83517 541.29 

THE MATRIX VALUES ABC, DE ETE. 

0.018182 0.042060 0.350738 0.032314 0.003923 


DODTRO = 2.32 DODTR3 = 190.90 DODTN = 328.70 
NUMBER OF ITERATIONS ARE:RIN 5 „FOR 0 12 

DODTN = 378.94 DMDT = 66.55 

RIN = 49.44 AIN = 110079.94 

VIN = 2423076.92 DVH20 = 1481890.78 

BEEN = 22.35 PDLTOT= 100.51 


THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS: 
630.00 623.42 8007.1 553.73 541.49 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.017470 0.041722 0.338296 0.032314 0.003923 


DODTRO = 2.11 DODTR3 = 219.69 DODIN = 378.94 
NUMBER OF ITERATIONS ARE:RIN 4 ,FORQ 12 

DODTN = 431.05 DMDT = 75.70 

RIN = 49.89 AIN = 5 

VIN = 2461538.46 DVH20 = 2272776.29 

PDL = 24.67 PDLTOT= 143.77 


THE TEMPERATURE PROFILE (TO,T2,T3,T4,T5) IS: 
640.00 632-27 614.93 555.02 541.69 

THE MATRIX VALUES AEB CDE ARES 

0.016849 0.041392 0.327034 0.032314 0.003923 


DODTRO = 1.89 DODTR3 = 249.70 DODTN = 431.05 
NUMBER OF ITERATIONS ARE:RIN 4 „FOR 0 12 

DODTN = 484.96 DMDT = 85.16 

RIN = 50.33 AIN = 111828.94 

VIN = 2500000.00 DVH20 = 4419910.85 

PDL = 26.95 PDLTOT= 257811 


THE TEMPERATURE PROEFILE (TO To T° sp n are: 
650.00 642.12 6 1 5۹57/7 541.90 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.016297 0.041070 0.316728 4+  ئى‎ ٣ 


DODTRO = 1.66 DODTR3 = 280.90 DODTN = 484.96 
NUMBER OF ITERATIONS ARE:RIN 5 ,FOR 0 12 

DODTN = 540.69 DMDT = 94.95 

RIN = 50.78 AIN = 112694.08 

VIN = 2538461.54 DVH20 = 8108853.38 

PDL = 29.21 PDLTOT= 448.55 


THE TEMPERATURE PROFILE (TO, DZ FE OR 
660.00 651.48 629.44 559.59 542.12 
IHE MATRIX VALUES A BCD E REE. 
. 0.015799 0.040754 0.307195 0.032314 0.003923 
DODTRO = 1.42 DQDTR3 = 313.30 DODTN = 540.69 
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APPENDIX A 


PROGRAM HOTSHIP 


The following FORTRAN program, as described in the text, 
is used to iterate the linearized heat transfer system of 
equations. Program inputs are U;,,, air temperature, and 
envelope emissivities. Several outputs, including the heat 
flux, q, and temperature profile values, are written to the 
screen and to files for use with the graphics program used 


æ compile figures 5 through 10. 


21 


10 


20 


100 


200 


THIS PROGRAM ITERATES A LINEAR HEAT TRANSFER MODEL OF AN AIRSHIP 
TO SOLVE FOR THE TEMPERATURE PROFILE AND THE RATE OF HEAT TRANSFER 


REAL MUI,MUMN,H1,H4,KMEAN,KINF,KSOV,KIN,KOUT,GRD,KFL,MUEUL,GREL 
REAL KRE , ۸٠ , ۱٢٢۷۸ ۰ E 11-7٧ 1ل له‎ ٢٢١ 


INPUTS ARE: THE ENVELOPE EMISSIVITIES, VELOCITY, AND AIR TEMP 
PI=3.1415927 

PRINT*,'DO YOU WANT TO INCLUDE THE RADIATION HEAT TRANSFER' 
PRINT*, "EFFECTS? (1=YES, 2=NO)' 

READ لا‎ 


>7 1 20 

PRINT*, 'INPUT THE EMISSIVITIES OF THE INNER AND OUTER' 
PRINT*, 'ENVELOPES, EIN AND EOUT' 

READ*, EIN,EOUT 


PRINT *,'INPUT U INFINITY (UI,KTS), AND T INF (DEG R)' 
READ *, UI,TI 


COMPUTSNGOIBESDIMENSIONS OPGPHECCZNHERSENVERBOPE 
WITH IHESE GIVEN OUTER ENVELOPE DIMENSIONS 
VTOT=4.0E6 

ROUT=66.5 

AOUT=1.4213E5 

[=399 23 


THIS LOOP INCREMENTS THE HELIUM TEMP IN 10 DEG STEPS 
DO 1200 K=1,14 

TO=TI+10.0*K 

M=1 

VIN=2.0E6*T0/520 

RINIT=50.0 
RTRY=(-3.0*VIN/(5.0*PI*(RINIT-3.0*L/5.0)))**0.5 
IF (ABS(RTRY-RINIT).LE.0.001) GO TO 200 
RINIT=RTRY 

IF (M.GE.20) GO TO 200 

M=M+1 

GO TO 6 

RIN=RTRY 


AIN= 2.0*PI*RIN**2+2.0*PI*RIN*(L-3.0*RIN)+PI*(RIN**2)*2.236068 
SPECIFYING THE PARAMETERS OF THE ENVELOPE MATERIALS 


KIN=0.03266/3600.0 
THIN=0.05/12.0 

KOUT=KIN 

THOUT=THIN 

IF UI=0, THEN FREE CONVECTION CONDITIONS EXIST 
IF (UI.EQ.0.0) GO TO 400 
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Figure 7 
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20 


Q vs Delta T 


V. CONCLUSIONS 


For this relatively simple concept design, heat loss is 
small enough to allow large changes in lift by heating the 
lifting gas with engine waste or exhaust heat or small 
liquid fuel furnaces. For envelopes with low emissivities 
produced as discussed, these results, and Appendix C 
calculations show that more than a 30% increase might be 
expected by use of engine waste heat alone under moderate 
power conditions. If steam injection is used, a 100% 
increase might be expected, with a significant time and fuel 
penalty required for evaporating the required amount of 
steam. 

By using heated gas, vectored thrust, dynamic lift 
and/or other lift augmentation systems, the airship as 
conceived might achieve more than a 100% increase in lift, 
which will allow a tremendous operational flexibility in 
cargo load capacity, unrefueled range and endurance and 
mission equipment carried. Obvious advantages over rotary 
wing aircraft in load carrying capacity include a much 
larger weight capacity, low fuel useage, longer range, and 
no rotor blast effects on takeoff or landing. Advantages 
over fixed wing aircraft are vertical landing and takeoff 


and lower fuel useage. 
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110.١ 5 


A. PROGRAM OUTPUT DATA 

Figures 5 through 10 show the program HOTSHIP results 
graphically for the parameters of interest. All results 
shown are for ambient air temperature of 60 deg F. 

Certainly this airship concept will have the same response 
to atmospheric temperature and density changes as a standard 
alrship. An additional minor effect will be the change in 
the inner envelope radius reguired to hold the same gross 
lift, which changes the insulative gap between the 
envelopes. 

It is quickly seen from Figures 5 through 7 that 
envelope emissivities will make a dramatic difference in 
heat loss. For organic materials, with emissivities greater 
than 0.9, the total heat loss rate is six or more times the 
convective heat loss alone (£ = 0.0). Even with envelope 
emissivity of 0.1, the radiative and convective heat losses 
are on the same order. 

The curve for the specific results Of s = 0.1 and 
Eout = 0.9 is shown, which is expected to correspond closely 
to the special case where the inner envelope only is vacuum 
metallized. Note that the heat loss for this case is less 


than if both envelopes have an emissivity of 0.2. 
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RIN = 51.22 AIN = 113552 4 
ئ٣‎ = 25 6923. 06 DVH2O = 35853817.74 
PDL = 30059 EDLIOT= 3 5 


fie TEMPERATURE PROFILE (TO,T2,T3,T4,T5) IS: 
670.00 614.45 587.96 3 6 566.08 
ØE MATRIX VALUES A,B,C,D,E ARE: 
0.008360 0.004045 0.000000 0.003231 0.003923 
DODTRO = 2:31. DODTR3 = 0.00 DODTN = 6647.57 


NUMBER OF ITERATIONS ARE:RIN 5 ,FOR Q 9 


DODTN = 7247.43 DMDT = 1 ے2‎ 

RIN: = 51.66 AIN = 114405.53 
وال 02310533162 = ج9‎ 20 5 KARRI 
PDL = 32.78 PDLTOT= 251905 


THE TEMPERATURE PROFILE (TO,T2,T3,T4,T5) IS: 
680.00 620.94 591.85 591.85 568.43 
THE MATRIX VALUES A,B,C,D,E ARE: 
0.008151 0.004014 0.000000 0.003231 0.003923 
DODTRO = 2.17 DODTR3 = 0.00 DODIN = 7247.43 
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THE TEMPERATURE PROFILE (TO,T2,T3,T4,T5) IS: 
630.00 8 618.84 541.81 540.00 

WHETMATRIX VALUES A,B,C,D,E ARE: 

0.015844 0.004172 0.314878 0.003231 0.000000 


DODTRO = 2.16 DODTR3 = 314.35 DODTN = 559.01 
NUMBER OF ITERATIONS ARE:RIN 4 „FOR O 11 

DODTN = 641.71 DMDT = 112.69 

RIN = 49.89 AIN = 5 

VIN = 2461533.46 DVH20 = 2272776.29 

PDL = 24.35 PDLTOT= 143.94 


TEMPERATURE PROFILE (T0,12,13,14,T5) IS:‏ كارن 
540.00 542.07 واو SA‏ 830225 640.00 

UHEOMATRIX VALUES A,B,C,D,E ARE: 

0.015248 0.004139 0.303883 0.003231 0.000000 


DODTRO = 1.95 DODTR3 = 360.30 DODTN = 641.71 
NUMBER OF ITERATIONS ARE:RIN 4 „FOR O 11 

DODTN = 728.44 DMDT = 127.92 

PIN = 50.33 AIN = 111828.94 

VIN = 2500000.00 DVH20 = 4419910.85 

FDL = 26.61 PDLTOT= 257.80 


BEESTEMPERATURE PROFILE (TO,T2,T3,14,T5) IS: 
650.00 629-31 636.31 542.35 540.00 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.014716 0.004107 0.293819 0.003231 0.000000 


DODTRO = 1.73 DODTR3 = 408.65 DODTN = 728.44 
NUMBER OF ITERATIONS ARE:RIN 5 „FOR O 11 

DODTN = 519-31" DMDT = 143.88 

RIN = 50.78 AIN = 112694.08 

VIN = 2538461.54 DVH20 = 8108853.38 

PDL = 28.84 PDLTOT= 450.24 


THE TEMPERATURE PROFILE (TO,T2,T3,T4,T5) IS: 
660.00 648.36 645.02 542.65 540.00 

UEESMATRIA VALUES A,B,C,D,E ARE: 

0.014236 0.004075 0.284508 0.003231 0.000000 


DODTRO = 1.50 DODTR3 = 459.49 DODTN = 819.31 
NUMBER OF ITERATIONS ARE:RIN 5 ,FORQ 11 

DODTN = 914.39 DMDT = 160.58 

RIN = 51.22 AIN = 113552.84 

VIN = 2576923.08 DVH20 = 35853817.74 

PDL = 31.04 PDLTOT= 1882.36 


THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS: 
670.00 657.40 653.70 542.95 540.00 
TEE MATRIX VALUES A,B,C,D,E ARE: 
0.013799 0.004045 0.275825 0.003231 0.000000 
DODTRO = 1.26 DODTR3 = 512.88 DODTN = 914.39 
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NUMBER OF ITERATIONS ARE:RIN 5 „FOR O 10 


DODTN = 267.75 DMDT = 47.02 
RIN = 47.61 AIN = 106500.01 
VIN = 2269230.77 DVH20 = 333175.90 
PDL = 12.55 PDLTOT= 30.61 


THE TEMPERATURE PROFILE (T0,T2 11٤ 
590700 584.87 592 540.87 540.00 

THE MATRIA VALUES A,B,€,D E ARES 

0.019363 0.004312 0.374875 0.093227 0.000000 


DODTRO = 2.92 DODTR3 = 153.44 DODTN = 267.75 
NUMBER OF ITERATIONS ARE:RIN 5 ,FOR 0 10 

DODIN z 334.67 DMDT = 58.77 

RIN = 48.07 AWN = 107205886 

VIN = 2307692.31 DV=70 = 453004.31 

PDL = 14.97 PDLTOT= . 89.38 


THE TEMPERATURE PROFILE )10,12,12 13 ٢٢ لا‎ 
600.00 593793 59. 524 541.08 540.00 
THE MATRIX VALUES ATB C D E ARER 


0.018230 0.004276 0.356389 0.003231 0.000000 


DODTRO = 2.74 DODTR3 = 190.35 DODTN = 334.67 
NUMBER OF ITERATIONS ARE:RIN 5 „FOR 0 10 

DODTN = 405.52 DMDT = 71.21 

RIN = 48.53 AIN = 108304.04 

VIN = 2346153.85 DVH20 = 655969.36 

PDL = 17.36 PDLTOT= 52.50 


THE TEMPERATURE PROFILE (TO, T2 ۳ ود‎ 7٤ 
610.00 603703 601.31 541:31 540.00 

THE MATRIX VALUES A,B,C,D, E ٣ 

0.017301 0.004241 0.340714 0.003231 0.000000 


DODTRO = 2.56 DODTR3 = 229.42 DODTN = 405.52 
NUMBER OF ITERATIONS ARE:RIN 5 „FOR Q 10 

DODTN = 480.27 DMDT = 84.34 

RIN = 48.98 AIN = £109195.42 

VIN = 2384615.38 DVH20 = 1071527.20 

PDL = 19,72 EDLIOT= 1227632 


THE TEMPERATURE PROFILE (TO T2 P Ta ۷ 
620700 5612 ٦1 610.09 ما‎ 540.00 

IHE MATRIX VALUES A,B,C D E ARET 

0.016516 0.004206 0.327048 0.003231 0.000000 


DODTRO = 2.36 DODTR3 = 270.72 DODTN = 480.27 
NUMBER OF ITERATIONS ARE:RIN 5 ,FOR 0 11 

DODTN z 559.01 DMDT = 98.17 

RIN = 49.44 AIN = 110079.94 

VIN = 2423076.92 DVH20 = 1481890.78 

PDL = 22.05 PDLTOT= 100.51 
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expansion of a fixed mass of Helium and the decrease in 
density of the air volume between the envelopes: 
me(lbm) = q(BTU)/20500(BTU/1bm) 


ehange 00101016 (1E) ~- (To - Tinf) 130,000 
TINE 


+ = = NE nh Cie hint) 


Mini 
“When steam injection is used to heat the inner volume, 
significant lift is produced by the volume the saturated 
steam displaces. Since the weight of the water is unchanged 


when it is vaporized, the change in lift produced by the 


steam is: 
home eM (steam eb) = change in VÔ 0۶ے‎ 
where een is the increase in steam volume from initial 


(Ting) to heated (T,) conditions. This volume is calculated 
for a specified temperature from: 

eae aa R/ [18.02 (SPVIDERT] 
using the ideal gas equation for Helium and steam table 
values from Subroutine Volume for the steam. 

The loss of air spacing between the envelopes due to the 
additional inner volume from the steam is not accounted for 
by the program. 

Output quantities are filed for use by the IBM graphics 


program DISSPLA. 
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Eg Fuel Flow vs Delta T (large scale) 


Figure 6: 


This concept has several significant disadvantages: 


- The double envelope will be heavier than conventional 
airship envelopes, although this factor is mitigated by 
not requiring additional ballonet envelopes. 


- A detailed design for the suspension and shape control 
of the envelopes has not been completed or tested. 


- Hot air ducting and/or a liquid fuel furnace will add 
` weight to the basic airship. 


? 


- A steam injection heating method requires a lot of water 
(14000 lb), energy to heat that water (7.32x10° BTU for 
a 50 deg F temperature rise = 360 lbm fuel), pays a 
weight penalty for the liquid film on the inner envelope 
(about 7000 lb), and requires some type of recovery 
system. 


- Avallability, toughness, and cost of a low emissivity 
fabric coating are unknown factors. 


- Some time is required for heatup or cool down during 
takeoff and landing operations where load changes are 
significant, depending on heating method used and heat 
source. 

Overall, the promise of operational flexibility and 
capability should lead to an expanded study of this concept, 
with the effort and financing required to accomplish large 
scale testing to prove the estimated and unknown design 


factors, and solve those engineering and conceptual problems 


which will arise. 
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1000 


1100 
1110 
1120 
21:0 
1140 


DODT=DODTN 

IF (UI.E0.0.0) GO TO 500 

GO TO 600 

CONTINUE 

EOVIVALENT FUEL BURN RATE AND LIFT CHANGE VALUES ARE 
DIRECTLY COMPUTED 


DMDT=DODTN*3600.0/20500.0 


DLIFT=130000.0*( (DELT/TI) )+( (2MEAN-TI)/ TI) =(VICI-Vin) “Rio: 


+*32.174 


VALUES BASED ON STEAM INJECTION HEATING ARE COMPUTED 
PDL=DLIFT*100.0/130000.0 


CALL VOLUME (TI,SPVI) 
VH201-VIN*T1*1545.0/(18.02*SPVI*14.7*144.0-1545.0*TI) 
CALL VOLUME (TO,SPVO) 
VH2OE=VIN*T0*1545.0/(18.02*SPVO*14.7*144.0-1545.0*TO) 
DVH20=VH20H-VH201 


DLH20=DVHA20*RH01*32.174-(VH20H*(TMEAN=TI)*REOMN*S2 0174 7۲ 
DLTOT=DLIFT+DLH20 
PDLTOT=(DLTOT*100.0)/130000.0 


PROGRAM OUTPUTS ARE PRINTED; FILES POR GRAPHICS USE ARE: WRI لتا‎ 


WRITE (*,1110) M,J 

WRITE (*,1120) DODTN,DMDT 

WRITE (*,1130) RIN, AIN 

WRITE (*,1140) VIN, DVEZO 

WRITE (*,1150) PDL, PDLTOT 

PRINT*, 'THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS:' 
WRITE (+. 1100) TO 12 73 ۶ 

PRINT*, ‘THE MATRIX VALUES A,B,C,D,E ARE: ' 

WRITE (*,1160) A,B,C,D,E 

WRITE (*,1180) DODTRO, DODTR3, DODTN 


WRITE (171100) VISTO, DELI DID PDL 50190010) 

WRITE(C2 1100) UT ,10,DELT DODTN FDL ٢٢٢ 

۷۶ ۲75 )3 7, 71 

WRITE(4,1100) TO,DELT,EIN,DODTN,DMDT,PDL 
FORHAT(F8.2,T10,F8.2,T20,F8.2,T30,F8.2,T40, F8.2,T50 7 F3.2 TOO EEE 
FORMAT(1X,'NUMBER OF ITERATIONS ARE:RIN',I2,' ,FOR Q ',I2) 


FORMAT(1X,'DODTN = ',F10.2,' DMDT = ',F10.2) 
FORMAT(1X,'RIN = ',F10.2,' AIN = ',F11.2) 
FORMAT(1X,'VIN = ',F10.2,' DVH20 = ',F11.2) 


DD 


SECTION 1: VALUES OBTAINED FROM THE BASIC PROGRAM 


INPUT THE EMISSIVITIES OF THE INNER AND OUTER 


ENVELOPES, EIN AND EOUT 
? 


.1,.9 

INPUT U INFINITY (UI,KTS), AND T INF (DEG R) 
2 

50,540 

NUMBER OF ITERATIONS ARE:RIN 5 ,FORQ 8 
DODTN = 41.16 DMDT = 7023 

RIN = 45.75 AIN = 102801.46 

VIN = 2115384.62 DVH20 = 25472.64 

PDL = 2.57 PDLTOT= 3.98 


MES TEMPERATURE PROFILE (TO0,T2,13,T4,T5) IS: 
350.00 548.80 548.61 540.29 540.16 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.031986 0.004468 0.549943 0.003231 0.003923 


DODTRO = 3.54 DODTR3 = 26.04 DODTN = 41.16 
NUMBER OF ITERATIONS ARE:RIN 5 ,FORQ 9 

DODTN = 89.71 DMDT = 15.75 

RIN = 46.22 AIN = 103737.93 

VIN = 2153846.15 DVH20 = 79283.17 

PDL = 5.14 PDLTOT= 9.50 


DE TEMPERATURE PROFILE (TO,T2,T3,T4,T5) IS: 
560.00 55776 557.36 S40.64 S40.35 

RHESMATRIA VALUES A,B,C,D,E ARE: 

0.025960 0.004427 0.471630 0.003231 0.003923 


DODTRO z 3.40 DODTR3 = 54.26 DODTN = 89.71 
NUMBER OF ITERATIONS ARE:RIN 5 „FOR O 9 

DODTN = 143.07 DMDT = 25.12 

RIN = 46.68 AIN = 104666.37 

VIN = 2192307.69 DVH20 = 121494.44 

PDL = 7.67 PDLTOT= 14.32 


IPS TEMPERATURE PROFILE (T0,T2,T3,14,T5) IS: 
57.0200 566 +0 995.17 2 540.56 

UNESMATRIX VALUES A,B,C,D,E ARE: 

BO 2913107004388. 0.428758 0.003231 0.003923 


DODTRO = 3.24 DODTR3 = 84.42 DODTN = 143.07 
HUMBER OH ITERATIONS ARE:RIN 5 ,FOR 0 10 

DODIN z 200.58 DMDT = 35.22 

RIN = 47.15 AIN = -‫ 19 

VIN = 2230769.23 DVHZO = 192264.72 

PDL _ = 10.18 PDLIOT= 20.63 


THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS: 
580.00 575.87 575.00 541.44 540.79 
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APPENDIX B 


PROGRAM SAMPLE OUTPUTS 


The following program outputs are used to derive the 


values quoted in the results and conclusions chapters. 


These are divided into the following sections: 


Section 


Section 


Section 


Section 


Section 


qu 
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5: 


Values obtained from the basic progran. 


Values obtained for T5 


To by setting 
to zero in the program. 


Values obtained for T, T, by, setting 이 
to zero in the program. 


Values obtained for T- = T by setting "E" 


3 
to zero in the program. 
Values obtained for increased envelope heat 
loss resistance by factors of 2 and 10. 


Units applicable to these data are: 


Distances/lengths in ft. 


Volumes in ft?; Areas in ft?. 


Heat flux in BTU/sec. 


Fuel flow rate in lbm/hr. 


Temperature in degrees R. 
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1150 
1760 
1180 


E 
1200 


۲۵٢۶۲۶۲۷1٦ CPDL ELO 2! PDLTOT= ',F11.2) 
BOPO -O ULL, FG-0,121,=6-o0,131,F9.6,141,F8.6) 


eri 1 20 للا‎ 57 232 DODIRI = ',F7.2,'" DODTN = !,F7.2/) 
THE PROGRAM IS RUN FOR OTHER INPUT VALUES IF DESIRED 

CONTINUE 

PRINT* 

BEJN 


PRINT*, 'DO YOU WANT TO ENTER OTHER VALUES OF UI,TI,EIN OR EOUT?' 
PRINT*, '(1=YES, 2=NO)! 

READ*, Z 

IF (Z.EQ.1.0) GO TO 10 

STOP 

END 

SUBROUTINE VOLUME (T,SPV) | 

THIS SUBROUTINE USES A SERIES OF LINEAR APPROXIMATIONS TO ' 
COMPUTE THE SPECIFIC VOLUME OF STEAM AT l ATMOSPHERE AT TEMP=T 


NEMO THEN 
SPV=2445.0+(T-500.0)*(867.7-2445.0)/30 


ELSEIF (T.LE.560.0) THEN 
SPV=867.7+(T-530.0)*(350.0-867.7)/30 


27 ۶ْ 590.0) THEN 
SPV=350.0+(T-560.0)*(157.17-350.0)/30 


۶ ۶۶۶ 9 0ئ‎ THEN 
Sev- lon Mai -590.0)Z(77.23-15717)/30 


ELSE CTE SSD) THEN 
5 20277 235 07-520 0405-00 


ELSE 
SPV=40.95+(T-650.0)*(23.15-40.95)/30 


ENDIF 


END 


56 


tabulated values at an intermediate temperature. A series 

of linear interpolations in Subroutine Volume is used to 

evaluate the specific volume of steam at a given temperature 

based on tabulated values from reference 12. Standard 

atmospheric pressure (14.7 psi) is used in all calculations. 
Since the linearized system of equations is of the 


simple form 


T = 


1.0 0 T, Ne و‎ 
l —1 O, OEE T4 0 
O 1-1 0 -C Tạ = =" 
0 0 1 =1 1 ič O 
HU UNO U q Tint 


a simple algebraic solution is used in the program for q, 
and then a substitution of new values is used to evaluate 
the temperature for the new iteration. 

Special computations are made if airship velocity or 
envelope emissivities are zero to avoid division by zero 
problems in the program. Results for airship velocities 
less than ten knots (but not zero) are less than accurate 
values due to the mixed convection conditions over the outer 
envelope. 

Equivalent fuel flow rate is based on 20500 BTU of 
useable heat energy per pound of fuel, which is obtainable 
from petroleum base fuels. Increase in lift due to heating 


Without steam injection is based on the sum from the volume 
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THE MATRIX VALUES A,B,C,D,E ARE: 
0.000000 0.004350 0.390713 0.003231 0.003923 


DODTRO = 3.01 DQDTR3 - 131.32 DODTN = 227.01 
NUMBER OF ITERATIONS ARE:RIN 5 ,FORQ 6 

DODTN = 296.14 DMDT = 52.01 

RIN = 47.61 AIN = 1 

VIN = 2269230.77 DVH20 = 333175.90  . 

PDL = 13.02 PDLTOT= 30.98 


Mae TEMPERATURE PROFILE (T0,T2,13,T4,TS) IS: 
590.00 0 2 542.12 541.16 

Ore MATRIX VALUES A,B,C,D.E ARE: 

92000000 0.004312 0.368715 0.003231 0.003923 


DODTRO = 2.82 DODTR3 = 169.74 DODTN = 296.14 
NUMBER OF ITERATIONS ARE:RIN 5 ,FORQ 6 

DODTN = 369.62 DMDT = . 64.91 

RIN = 48.07 AIN = 3 

VIN = 2307692.31 DVH20 = 453004.31 

SDL = 15.52 PDLTOT= 39.76 


EM TEMPERATURE PROFILE (IO,T2,T3,I4,T5) IS: 
600.00 0 598.42 542.64 541.45 


THE MATRIX VALUES A,B,C,D,E ARE: 
0.000000 0.004276 0.350756 0.003231 0.003923 


DODTRO = 2.62 DODTR3 = 210.61 DODTN = 369.62 
NUMBER OF ITERATIONS ARE:RIN 5 ,FORQ 6 

DODTN = 447.46 DMDT = 78.58 

RIN = 48.53 AIN = 108304.04 

VIN = 2346153.85 DVH20 = 655969.36 

PDL = 17.98 PDLTOT= 52.85 


DEESTEMPERATURE PROFILE (TO,I2,T3,I14,T5) IS: 
610.00 610.00 608.10 543.20 541.76 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.000000 0.004241 0.335507 0.003231 0.003923 


DODTRO = 2.41 DODTR3 = 254.02 DODTN = 447.46 
NUMBER OF ITERATIONS ARE:RIN 5 ,FORQ 6 

DODTN = 529.68 DMDT = 93.02 

RIN = 48.98 AIN = 2ه‎ 

VIN = 2384615.38 DVH20 = 1071527.20 

PDL = 20.41 PDLTOT= 77.02 


THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS: 
620.00 620.00 62 543.79 542.08 
AS MATRIZ VALUES A,B,C,D,E ARE: 
0.000000 0.004206 0.322198 0.003231 0.003923 
DODTRO = 219 1013-5 3009.06 DODTN = 529.68 
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THE MATRIK VALUES ABC DE REET. 
0.012216 0.004350 0.000000 0.003237 07002773 


DODTRO = 3.37  DODTR3 z 0.00 DODTN = 1688.06 
NUMBER OF ITERATIONS ARE:RIN 5 ,FORQ 8 

DODTN = 2188.42 DMDT = 384.31 

RIN = 47.61 AIN = 106500.01 

VIN = 2269230.77 DVH20 = 333175.90 

PDL = 12.02 PDLTOT= 30.21 


THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS: 
590.00 565.09 555.66 555.66 548.58 
THE MATRIX VALUES A,B,C,D,E ARE: 
0.011398 0.004312 0.000000 0.003231 0.003923 
DODTRO = 3.27 DODTR3 = 0-00 DODTN = 2188.42 


NUMBER OF ITERATIONS ARE:RIN 5 EOS 


DODTN = 2705.01 DMDT = 475.03 
RIN = 48.07 AIN = 107405.63 
VIN = 2307692.31 DVH20 = 453004.31 
PDL = 14.42 PDLTOT= 39.00 


THE TEMPERATURE PROFILE (T0,T2,T3 127 TS) 200; 
600.00 570732 559255 22935 550.6 


THE MATRIX VALUES A,B,C,D,E ARE: 
0.010763 0.004276 0.000000 0.003231 0.003923 
DODTRO = 3.17 DODTR3 = 0.00 DODTN = 2705.01 


NUMBER OF ITERATIONS ARE RIN 5 ,FORO 8 


DODTN = 3235.65 F DIDI = 5060.21 

RIN = 48.53 AIN = 108304.04 
VIN = 2346153.85 DVHZO = 655969.36 
PDL = 16:21 6U i0 32.20 


THE TEMPERATURE PROFILE (TO,T2,T3,T4,T5) IS: 
610.00 576.87 563.15 563.15 552.69 
THE MATRIX VALUES A,B,C,D,E ARE: 
0.010249 0.004241 0.000000 0.003231 0.003923 
DODTRO = 3.07 DODTR3 = 0.00 DQDTN = 3235.65 


NUMBER OF ITERATLONS +۹۹۷۹٦ 


DODTN = 3778.79 DMDT = 663.59 

RIN = 48.98 AIN = 109195.42 
VIN = 2384615.38 DVH20 = 1071527.20 
PDL = 19.17 PDLTOT= 76.66 


THE TEMPERATURE PROFILE (TO,T2,T3,T4,T5) IS: 
620.00 582.093 567.03 567.08 554.82 

THE MATRIX VALUES TA P CAD P RE: 

0.009818 0.004206 0.000000 0.003231 70 0909923 


DODTRO = 2.95 DODTR3 = 0.00 DODIN = 3773279 
NUMBER OF ITERATIONS ARE:RIN 5 „FOR Q 9 

DODTN = 4333.20 DHDT = 760595 

RIN = 49.44 AIN = 110079 . 4 
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The radiation equations are basically analytical, and 
should be accurate to within the uncertainties in 
emissivities and temperatures. 

Overall unknown factors include additional conductive 
and convective losses in the lower portion of the airship 
due to the anchoring of the inner envelope, and losses in 
ducting and structural components. 

Although the combination of these uncertainty factors 
is: 

Go (1.4)(125)(1,25) = 3.15 
the major factor is the convection between the envelopes. 
It is believed that an airship constructed as conceived will 
have convective heat losses well within a factor of two 
times these results. Radiative heat losses, restricted to 
the assumptions stated, will be within a few percent of 


these results. 
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Figure 5 


When this linear system is solved, these new solutions 
are resubstituted into the equations to find the heat 
~ ster coefficients and radiative heat transfer values for 
the new iteration. This entire problem is iterated until 
convergence for g is obtained. 

Since several of the "constants" in the linearized 
system depend on the temperatures (Rin: hi, h,, and ا‎ and 
particularly since the radiative heat transfer terms are 
fourth order in temperature, convergence of this iteration 
is expected to be a major problem. A close initial guess of 
the temperature values, however, produces convergence up to 
critical envelope emissivity values of more than 0.9. This 


is considered an adequate range for emissivities, and an 


adequate solution for this model. 


E. THE COMPUTER PROGRAM 

Program HOTSHIP utilizes these equations and completes 
the iterations with a convergence criterion of q within 0.1 
BTU/sec. From input values of airship velocity, Tinf» To 
and envelope emissivities, the temperature profile and heat 
loss rate are computed. From these values, the increase in 
lift and equivalent fuel burn rate are directly computed. 

Properties of air and Helium with little temperature 
dependence are evaluated by a linear interpolation from 
tabulated values obtained from reference 8. The ideal gas 


equation is used in density change computations based on 
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NUMBER OF ITERATIONS 25:11 ۶ EOE OT 


DODTN = 545.87 DMDT = 25.36 

RIN = 49.44 AIN = 110079 . 4 
VIN = 2423076.92 DVH20 = 1481890.78 
PDL = 22.24 PDLTOT= 100.71 


THE TEMPERATURE PROFILE (TO,T2,T3,T4,T5) IS: 
630500 62 6 619.06 943.91 542.14 

THE MATRIX VALUES A,B € ETARE: 

0.015939 0.0041727 07376473 70.003237  7>٭‎ 


DODTRO = 2.16 DODTR3 = 308.42 DODTN = 545.87 
NUMBER OF ITERATIONS ARE:RIN 4 ,FOR Q 11 

DODTN = 626.30 DMDT = 109.98 

RIN = 49.89 AIN = 5 

VIN = 2461538.46 DVH20 = 2272776.29 

PDL = 24.55 PDLTOT= 143.83 


THE TEMPERATURE PROFILE (TO,T2,T3,T4,T5) IS: 
640.00 630.42 627.83 544.48 542.46 
THE MATRIX VALUES A,B,C,D,E ARE: 
0.015341 0.004139 0.305484 0.003231 0.003923 
DODTRO = 1.95 DODTR3 = 353.46 DODTN = 626.30 


NUMBER OF ITERATIONS ARE:RIN 4 ,FOR Q II 


DODTN = FLOR 1) = 124.79 

RIN = 50.33 AIN = 111828.94 
VIN = 2500000.00 DVH20 = 4419910.85 
PDL = 26.83 PDLTOT= 25739 


THE TEMPERATURE PROFILE (TO,T2 1 159). 
650.00 639250 66 5 545.08 542.79 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.014808 0.004107 0.295422 0.003231 02003923 


DODTRO = 1.73 DODTR3 = 400.86 DODTN = 710.59 


NUMBER OF ITERATIONS ARE:RIN 5 „FOR O ll 


DODTN = 198.84 DMDT = 140.28 

RIN = 50.78 AIN = 112694.08 
VIN = 2538461.54 DVH20 = 8108853.38 
PDL = 29.08 PDLTOT= . 449.13 


THE TEMPERATURE PROFILE (TO,T2,T3,T4,T5) IS: 
660.00 648.58 645.32 545.72 543.13 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.014327 0.004075 0.286114 0.003231 0.003923 


DODTRO = 1.50 DODTR3 = 450.71 DODIN = 798.84 
NUMBER OF ITERATIONS ARE:RIN 5 ,FOR Q 11 

DODTN = 891.11 DMDT = 156.49 

RIN = 51.22 AIN = 113552.84 

VIN = 2576923.08 DVH20 = 35853817.74 


PDL 31.320 5 1876-12 
THE TEMPERATURE PROFILE (T0,T2,73,14,75) Is: 
670.00 657.64 654.04 546.38 543.50 
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NUMBER OF ITERATIONS ARE:RIN 5 „FOR O 12 


DODTN = 598.22 DMDT = 105.05 
Di = 51.22 RIN = 113552.84 
VIN = 2576923.08 DVH20 = 35853817.74 
DL = DOT 1872.96 


MARE TEMPERATURE PROFILE (10,T2,T3,14,T9) IS: 
670.00 660.84 636.64 561.68 542.35 

BEER MATRIX VALUES A,B,C,D,E ARE: 

0.015348 0.040446 0.298304 0.032314 0.003923 


DODTRO = 1.17 DỌDTR3 = 346.92 DODTN = 598.22 
NUMBER OF ITERATIONS ARE:RIN 5 ,FOR 0 13 

DODTN = 657.62 DMDT = 115.48 

RIN = 51.66 AIN = 53 

VIN = 2615384.62 DVH20 = 大 大 大 大 大 大火 六 大 大 六 

TU = 33.62 PDLTOT= : | -807.37 


THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS: 
680.00 670.19 643.79 563.83 542.58 
THE MATRIX VALUES A,B,C,D,E ARE: 
0.014937 0.040145 0.289940 0.032314 0.003923 
DODTRO = 0.91 DODTR3 = 381.83 DODTN = 657.62 
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APPENDIX C 


COMPUTATION OF AVAILABLE ENGINE HEAT 


From publications on powerplant characteristics, such as 
the engine manufacturer's specifications or Aviation Week, 
some basic trends for turboshaft and turboprop engines such 
as would be used for propulsion on this airship concept are: 

exhaust gas temperature = 600 deg C = 1100 deg F 
mass flow rate (lbm/sec) = shp/110 (lbm/sec) 

It is estimated that two 1000 shp engines (at least) 
would be used on the concept design. If the exhaust 
temperature could be reduced by 500 deg F in a heat exchange 
system without significantly effecting engine performance, 


then the exhaust heat available from each engine is: 


G (mass flow rate) Cp4i, (500) 


(1000/110) (7.7/32.174) (500) BTU/sec 


1088 BTU/sec 

Certainly almost twice as much heat (a 1000 deg F drop 
in temperature) is available if the exhaust is ducted 
directly into an air volume. 

Some additional heat is available from cooling the 
engine casing, but an estimate of quantity is not completed. 

Approximately one third of these heat values should be 
available at low or idle power, and about one half is 


expected to be available at moderate power. 
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APPENDIX D 


E CALCULATION OF VALUES FOR STEAM INJECTION 


1. Estimation of Water Film Weight 
From Gebhart (reference 10, pp. 435-436) the average 
thickness of a liquid water film on a large vertical plate 
will be about 0.005 inches. Based on this thickness, 
WelghtH2o film = Area (thickness) P g 


1.5x10>ft“ (.005 ft/12) (60.81b/ft?) 


3800 LDE 

From the additional film thickness at the top and bottom 
of the envelope, the total water film weight is expected to 
double this value, or as much as 7000 lbf. 
2. Mass of Steam Evaporated 

At 110 deg F, the mass of water vapor in 2x109ft? of 
saturated Helium is: 

Moo V/Viao = 2x10% /265.1 = 7544 lbm 
3. Energy to Evaporate Steam 

For the mass of water evaporated at 110 deg F, 

heat required = (970BTU/lbm)(7544lbm) 

- 7.32 x 109*"BTU 

or in terms of equivalent fuel (kerosene), 
= 7.32x10° BTU/(20500 BTU/lbm) 


"kerosene 


360 lbm 


If engine heat of 2000 BTU/sec is used, then the time 
required to evaporate this much steam is: 
time = [7.32x109/2000](1 hr/3600 sec) = 1 hr 
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